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ABSTRACT 


Computer Simulation can be of considerable assistance 
inggcheseinclusion oftlispace “aspmag* variable inkecological 
systems. in Chapter 2, the major concepts of ecology and 
ecological modeling are defined. Simulation is then 
ennsraered fas a* etool Gfore’studying secological g#nodeliss 
Chapter 3 reviews in detail spatial dynamics in community 
ecology. Examples from the literature are given. 

The fourth chapter describes ECOSIM, which allows a 
user to model the behavior of an ecological community in 
both time and space. The distribution, mobility, breeding, 
predation, and hunger characteristics of each species can be 
specified. Selective barriers to mobility can be 
established. Model definition, simulation monitoring, and 
editing facilities are provided. Compact but extendable 
data structures are employed. Output capabilities include 
graphical display of the simulation space and the optionvof 
plotting the population of given species against time. 
States of the model can be stored and later retrieved for a 
new Simulation run. 

In the final chapter, a simulation study illustrating 
some features of the system is presented, the system is 
evaluated, and some improvements to it are suggested. An 


Appendix contains a User's Manual. 
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CHAPTER 1 


INTRODUCTION 


The methodologies of differential equations and 
statistics have been widely used to study population growth 
in ecological systems. Multi-species or community growth 
has also been studied uSing statistical methods, but without 
much success because of the proliferation of relevant 
parameters’ (Bartlett, 1973).* As May (1976) obsérves;,;;one of 
the most efficient ways to study multi-species environments 
is to build a model including parameter values based on 
observation or intuition and then use a computer to simulate 
the dynamics of the system. 

Work has also been done on the patterns of organisms at 
one moment in/time (Pielou, 1969) but so far little progress 
has been made, even for single species populations, towards 
combining the study of how populations vary in time and in 
space. By modeling a two-dimensional environment, computer 
Simulation can be of considerable assistance as the 
inclusion of space as a variable in ecological systems. 

FCOSIM is an interactive programming system for use in 
the simulation of Spatial dynamics in community ecology. It 
allows a user to model the behavior of a community of 
species in both time and “space. The environment is 
represented aS an unbounded two-dimensional grid of 
hexagonal cells. The user can specify predation, mobility, 


breeding, and hunger characteristics for any number of 
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species. 

In the next chapter, the major concepts of ecology and 
ecological modeling are defined. Simulation is then 
constdencdgeas PalGtool for studying “ecological models. 
Chapter 3 reviews in detail spatial dynamics in community 
ecology. Dispersion, aggregation and Space-time 
interactions are discussed. Examples from the literature 
are given. 

ECOSIM is the subject of Chapter 4. The system is 
briefly described and the characteristics of the main 
modules are defined. The required data structures are 
presented in detail. There follows a definition of the main 
variables included in the model, and their use. The system 
requirements are then described and an evaluation of its 
performance iS given. Chapter 5 describes how ECOSIM is 
used to simulate a real-life problem, the system is 
evaluated and improvements to it are suggested. A User's 
Manual is included in the Appendix to provide detailed 


information on the use of the system. 
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GRAPTER TZ 


MODELING AND SIMULATION IN ECOLOGY 


2-1 ECOLOGICAL MODELS 


2e 121 Introduction 

From the Greek word otkos ("house"), ecology is the 
biological discipline studying the relationship between 
organisms and environment (Sacchi §& Testard, 1994). The 
problem of distribution and abundance of organisms is also 
of great interest to biologists (Krebs, 1972). the 
environment comprises the physical, chemical, and biological 
characteristics of the surroundings. 

Populations have been studied since the 19th century. 
In 1798, Malthus published the first book on demography. He 
studied the relationship between increase in population of 
organisms and increase of food supply. Some of the 
questions he raised remain unanswered today. 

Attempts to define mathematical relationships between 
organisms and environment have been partly successful. The 
ebassretiworky.ist? ithatecofenlLotkatecand: fYolterratein 1926 
(VYolterraye1i937) 6 pwAlthoughwmanysecologists have criticized 
their theory as simplistic and inadequate, the approach is 
Stillpiwadely guseds Inclusion of stochastic variates has 
contributed to the realism of more recent formalizations. 
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populations in an environment, a community. 


221.2 Population Ecology 

A population can be defined as a group of organisms of 
the same speciesteoccupying® astparticular space at a 
Pale culagytime e(hruebs ae 197 Ayes foretaseei so tatcaseunrtsthot 
interbreeding individuals (Emlen, 1973). Population ecology 
has been extensively studied since the work of Lotka and 
Volterra. Methods to estimate population growth, £er 
example, or prey-predator interactions are numerous. The 
basic parameters of a population include density, birth 
rate, death rate, immigration, emigration, and age 
distribution. With a variety of methods one can estimate 
the present or past values of these parameters. But the 
main purpose of the analysis is to predict phenomena such as 
the growth of a population or the effects of insertion of a 


predator in a species’ environment. 


221.221 Population Growth 

The growth of a population can be analyzed 
mathematically. This section treats different techniques, 
the parameters involved and the way to specify values for 
the different variables. 
If we first assume that 
there is no overlap between successive generations, so that 
time. can be expressed asa discrete variable (such cases 


exist in nature), we can express a population growth as: 
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Nee) Sone [1 eaer (Aenea ye] 
where 
Nee ~epopulation sizesateqenerationsteds: 
Ny = population size at generation t; 
r =ehnetqreproductive rate; 
K = "Carrying capacity of the environment", the 


population density that the environment can 
support. 
The population at time t+1 will thus be the population at 
time t (N;), Increased or decreased by a factor determined 
by deaths and births (rN,), and decreased by a factor which 
is a function of the density of the population (CN{/K). 
Continuous deterministic gqrowth. Generations usually 
overlap, making differential equations more suitable for 
representing population growth. We must first make the 
following assumptions: 
(1) breeding occurs uniformly throughout the year with 
no peaks; 
(2) breeding is independent of sex and genotype; 
(3) «the “influence “of density on»the rate of /increase 
operates without time lag. 
Population growth can then be expressed as; 
adN/dt = rN (1 - N/K) 
The solution of this differential equation is 
arb 
Ny = Ka7oGCiatec ) 
where a is a constant of integration defining the position 


of the curve relative to the origin. This type of equation 
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hasebecn® used/sto! modely tthe lgrowth*fofs' sone? Ycontrolled 
laboratory populations, yeast and protozoa (Gause, 1934) for 
exanple, Stand’ isome’sother® speciesrilikes \beetlesvand Eruit 
flies. But most natural populations do not seem to grow 
according to such eguations. 

Attempts were made to reduce the simplicity of the 
assumptions. Assuming that a change in WN _ simultaneously 
Gelcetse "pb" and.ed 91s notirealistic. “In ceality, time wags 
occur. For example, a high population level at time t may 
produce fae shortage ~of “food "at*tine ttt “and “ardiminution® rn 
thesbirth “patestat wrzs Hassel (1975) proposed a new 
discrete model for population growth which takes into 
consideration the big density-dependence which as 
Significant only with high population density. 

It was assumed for analysis that breeding occurs 
continuously throughout the year with no. peak. ieee acta, 
most animals have seasonal reproductive peaks so that r is 
not constant. Death and birth rates also vary by sex and 
age. The growth of a population is thus influenced by its 
sex and age structure. To include these factors and to _ be 
able to consider different age groups, two approaches are 
used: the life table and the Leslie matrix. 

An example of a life table is given in Fig. 2.1. This 
table describes the mortality schedule of a population. The 
tetatic! alife “tablieraisical culated fontthetbasrsrfoteatcross 
section ofladpopulationtat avspecifici. tine; whilega “cohort! 


life tablevis’ calculated®.on ‘the basis’ of a’ scohort (of 
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organisms followed throughout life. By calculating the net 
Eeproductive bated (Rk),  which=\is, the»! expected® numbertbof 
female offspring to which a newborn female will give birth 
in her lifetime, and the mean generation time (T), which is 
the average age at which females produce offspring, we can 
determine the population growth: the population WEL 


inerease) byjamfiact orsofcRtevery Tiyears. 
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0 142 1000 BS eee) Se ryan) 
1 62 437 198 0.453 1596 
2 34 239 98 0.410 Zevhd 
2 20 141 a2 Oa227 2.34 
4 (15.5) 12 109 a2 0.294 1.88 
11 77 31 0.403 124d 
6 (6.5) 1 46 32 0.696 CPP Ee 
A Z 14 0 0.000 eoO 
8 2 14 14 1.000 0.50 
9 0 0) 228 “8 “Saree aoe = 


(1) Age (year) x; 

(2) Observed number barnacles alive each year; 

(3) Number surviving at start of age interval x, ly; 

(4) Number dying within age interval (x,x+1), d,;3 

(5) Rate of mortality, qx; 

(6) Mean expectation of further life for animals alive at 
Start of agerxy @x5$ 

1 Estimated number alive. 


mm em me ee me me ee we we we we ee we we ee we ww ww we wwe wa ww ae ae iw ae oe ee es we = 


The Leslie method uses matrices to introduce the 
Variable’ages The structure of: the model representation 
allows calculation of the number of individuals of each age 


group at every time step. Also considered is the lag 
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betweentather tine of birthvofeanvindividualtand thertime he 
enters#tthe breedingyrpopulation: This method is thus 
convenient when a life table is available; the results are 
sometimes ciosely related to laboratory populations. 
Stochastic growth. The models discussed so far were 
deterministic. Biological systems are more realistically 
treated in probabilistic terms. We can often work with the 
probability that a female will have a litter or that a 
predator will kill a prey... A simple model (Pielou, 1969), 
assumes that if 
B(N) = birth rate and D(N) = death rate, then 
B(N) must decrease as N increases, 
B (Nanay —-ben cand 
D(N) must increase as N increases, 
D(N}eeea, +b, N. 


SO, —LOmSDELrth; 


P(N --> N+1) oc NB(N) = a,N - b,N2, 
and for death, 

P(N --> N-1) cc ND(N) = a,N + b,N2. 
Thus 

dN /dte=t (a, sa, )Ne- (bpetabsy Ne: 


Even though stochastic models are ‘more «realistic, the 
introduction of more variables makes the analysis very 
ditfiacuit. 

Most of the models described above predicted life 
cycles of some populations of simple organisms. But factors 


like large numbers of species, the uniqueness of each 
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individual and its continual adaptation to the environment, 
immigration and emigration, territorial behavior, and 


climatic fluctuations were not considered. 


22142.2 Species Interaction 

There are three types of interaction between two 
species. Predator-Prey relationships (+,-) occur when the 
growth rate of the predator species is increased by the 
occurence of the prey species and inversely. If the lLotka- 
Volterra equations are used many unrealistic assumptions 
must be made. Age and sex differences are again ignored. 
The effects of interactions bet ween species are 
instantaneous. The system is closed. Despite these 
restrictions, good results have been obtained, but mostly 
for laboratory populations. 

Competition (-,-) occurs" when thesgrowth rates» of» the 
two populations are depressed by the co-occurence of the two 
species. If the animals belong to the same species the 
eanpecitioneris! intraspecific ye interspecific competition 
exists between different species. Equations describing this 
phenomenon fic only a few laboratory populations. 
Transferring the results of laboratory work to field 
populations has proven difficult. 

When both growth rates are increased by the cooccurence 
Oba thedetwopsspéecies,. thes nellationshi pefispnutualasticaor 
commensal (+,+). Commensalism is "the relationship between 


two organisms when one is benefited" (Emlen, 1973, p.319). 
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The equations for competition remain valid for commensalism 
and mutualism. These relationships can be treated as 


negative competition. 


2.1.2.3 Population Regulation 

In reality, the growth potential of a species is never 
reached. Ecologists thus believe in some form of regulation 
to allow for more food and more shelters and less attraction 
for predators. This mechanism is not considered voluntary 
but a consequence of two types of limiting factors: density 
dependent factors, which are inversely related to population 
size, and density independent factors. Changes in the 
environment, availability of resources, predation, space and 
Spacmng jand “social factors like crowding all (contributes to 
population regulation and are not independent of each other. 
There has not yet been any attempt to introduce parameters 
to account for all of these factors ina systen. 

Fluctuation frequently occurs in population growth. 
Although fertility and survival rates are the main causes, 
delayed density dependent response and predator-prey 


interaction also influence the stability of a population. 


Zod 3aCommunity’ Ecology 

It is difficult to adequately represent an environment 
for a given population. The huge number of parameters and 
the fact that the populations coexist with many other 


species complicate the process. WThe. «Ste pe sup stom thres 
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dimensions introduces not only a confusing proliferation of 
parameters, but also a qualitative change in the dynamical 


conpkextey "e(MaypedoeG, ps69)2 


2.1.3.1 Mathematical Representation 


A community is an association of ranteracting 
individuals living together in the same environment. 
Mathematical analyses of population growth within a 


community are rare. The system is usually represented by a 
generalization of the Lotka-Volterra equations. For 
example, a community of competing species could be described 
as (Pianka, 1976) 


aNivat ={rpNoe (KysoNg=p oy Ny) /K 


where 
a = number of species; 
Egy = catexgor increase of species i; 
Ky, = carrying capacity of species i; 
Ni = population density of species i; 
Xj} = competition coefficient. 


Although the Lotka-Volterra equations have many flaws, 
they do provide a conceptual framework. Progress has been 
made towards defining population growth in a community. But 
most analyses of such systems are: still qualitative rather 


than quantitative. 
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2.1.3.2 Species Diversity and Abundance 

One of the most important properties of a community is 
the numberiotatspecitess it scontains. This number should 
contribute to a measure of the "diversity" of the community. 
But if one species present in the community is much more 
numerous than another, this disparity should be reflected in 
audiversitytindex. @fAn example ofsesuch tan bindexigiseethe 
Shannon-Wiener "entropy" measure: 

H = a p, log(p,) units, 

Wiere @#p4e*is 'othe peapercies of the i'fth species in the 
community and s is the number of species. p; can be the 
proportion of trees or animals or, in the case of plants, 
the relative weight of each species. 

Most ecologists agree that there exists a relationship 
between the diversity of a community and its stability. In 
theory, a diverse community “ought" to be stable. If there 
are many prey for a single predator, the predator population 
will less likely go extinct, and no prey population will 
grow indefinitely since there are many predator species. 

But it is not yet known how many species can coexist in 
a given environment. Two processes are involved. The 
number of species tends to decrease with competitive 
exclusionise)Gause” (1934) hassobserved (that pYliicass aniresult 
of competition two similar species ‘scarcely ever occupy 
similar niches, but displace each other «in such ©a manner 


that each takes possession of certain peculiar kinds of food 


andeemodes® of) (life sine whichvitshas*an advantage over its 
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competitor", On the other hand, the number of species tends 
to increase with immigration and evolution. 

A dominant species is one which is highly successful 
ecologically and has a great influence on the growth of the 
other species (Krebs, 1972). «Dominance is thus «related to 
the concept of species diversity. Although dominant species 
in a community are usually assumed to be ecologically 
constant, it happens frequently that a dominant species is 
largely a matter of chance events. Hence no mathematical 


model exists to cope with dominance. 


2.2 SIMULATION 


2.2.1 Introduction 

It is often difficult to analyze complex dynamic 
systems by means of mathematics. Equations describing 
changes in behavior may be intractable, or mathematical 
formulations may not exist at all. The complexity of large- 
scale, real-life Situations contributes to the 
oversimplification or lack of understanding responsible for 
most failures in attempting to mathematically define natural 
systems. 

Simulation can be used to analyze problems which are 
too complex to be defined or solved in terms of mathematics 
only. There are dozens of books on simulation and each of 
them has its own definition of the subject. Merging the 


concepts together, the following was obtained. 
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A system is a set of entities united by some form of 
interaction or interdependence (Mize, 1968). There are 
natural systems and artificial systems. A system is open if 
it can exist only in a particular environment and adaptive 
if it reacts to changes in environment. A model is an 
abstraction or a representation of a real world system. A 
simulation is the process of reproducing the system behavior 
using the model. A model can then be a system of equations 
and a simulation would be used to solve that system on a 
computer. This is referred to as a "mathematical" model, or 
as; Hall. andgdaye (1977) caldtmitjatan SNanalyticalve model: 
While modelling establishes the relationship between a real 
system and its model, simulation deals mainly with the 
relationship between the computer and the model (Zeigler, 


1976)%3 


2.2.2 Computer Simulation 

The first step in studying a system by computer is to 
build a model. This@Pinegdel @is #speciited@eby Barseteot 
variables representing the principal features of the real- 
life system and a set of computer instructions representing 
the laws or decision rules that determine how these features 
are modified as time progresses" (Smith, 1968). The systems 
involved may be of business, economic, social, biological, 
physical, or chemical type. A model is replicatively valid 
if it adequately reproduces the behavioral data from _ the 


real world system. A predictively valid model is one which 
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can accurately predict the future behavior of the real 
system. The use of simulation can make the model easier to 
understand, and can provide hints as to the causes and 
effects of the process. Simulation usually speeds up the 
analysis, provides a basis for system modifications, and 
facilitates experimental manipulation. Guides can be 
provided for predictions and forecasting. Control is 
possible over the time scale and sources of variation. The 
model will then enable the modeller to learn how change in 
State occurs, and to predict and control changes. 

R@Pnodel can be. deterministic sor stochastic,» Tti4as 
deterministic when all entities bear fixed mathematical or 
logical relationship to each other (Fishman, 1973). When 
the future is influenced by factors whose effect cannot be 
foreseen exactly, the model is stochastic. We can then 
describe the relationships between cause and effect with 
probabilities. 

Rienode! eisestaticeiffitedees tot involvesitiness sitnis 
dynamic when time is included as a variable. There are 
continuous time models (such as those based on differential 
equations), discrete-event models and discrete-time models 
(Zeigher?,eri1976) . A continuous time model is used when 
changes in the states of the model occur continuously as 
time evolves. in *avediscretestimeounodel ,Staretclocktais 
simulated and used to count time. In a discrete-event 
model, even though time flows continuously, state changes 


Can occur only in discontinuous jumps. A ymodel «is, tine 
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varying if time enters explicitly as an argument of the 
rules of interaction. Otherwise, the model is time 


invariant. 


2.2.3 Programming Considerations 

The construction of a computer simulation involves the 
following steps (Emshoff and Sisson, 1970; Ziegler, 1976): 
1. Define the problem. Study the real system to see which 


factors affect it; find the appropriate time scale. 


Zs) Analyze the data requirements and availability of 
sources. 
3. Define the experimental frame, which "characterizes a 


limited set of circumstances under which the real system 
is to be observed or experimented with" (Ziegler, 1976, 
pe 30)./ A subset of input-output behavior of the real 
system is associated with the experimental frame. 

LAY Defanetethe “base® tnodelyoliwhich*¥isievalidtiin, ali the 
allowable experimental frames. The base model accounts 
for all the’ input-output behavior of the real system. 

Beenpecreate ea “lumped modéely;weeatsubsysten of ftheshbase snodel 
which corresponds to the experimental frame of current 
interest. 

6. Validate the model. A common validation method is to 
compare results from the simulation to historical data. 
Because the usual purpose of a model is to predict the 
fuguremmeresultse ofimthetdpiferent simulations must, at 


some point, be checked with what actually happened. 
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7. Design simulation experiments. Specific combinations of 
manipulatable variables are used in the model to mimic 
real observations of the systen. 

8. Analyze the results. The simulation runs are evaluated 
and weaknesses are corrected if possible. 

Programming a simulation touches many aspects of 
computing science, including creating, filing, and 
destroying records, searching lists, and generating random 
variates. Classes of objects within a system must be 
defined, adjustments performed as conditions within the 
system vary, and objects related to one another and to their 
common environment. 

Consideration has to be given to random number 
generators. Random variables are used to represent factors 
which Eluctmate* (in) “an Sunpredictable™ buts statistically 
describable way. One can use a generator to simulate almost 
any discrete or continuous probability distribution. 

Another concern is whether the model is interactive or 
not. An interactive model is one that interacts with its 
user. Some of the reguirements far an interactive 
Simulation package include (Sohnle, Tartar, & Sampson, 
sagspyea*tiexvbleyvconmmand Language; easy \erore thee wsenys to 
understand; graphics and output capabilities so that the 
user can visualize the model; user capability to specify the 
length of the simulation, to save and restore the results 
and to restart the simulation; provisions to modify the 


model characteristics and parameters; and painless error 
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recovery. 

The language used for the actual coding of the 
Simulation program needs good data structure capabilities, 
conditional and unconditional branch and iteration 
Capability, and of course logical and arithmetic operations. 

Engineering design, econometrics, nuclear reactor 
design, military Operational research, and management 
science are all sources of models suitable for simulation. 
Many management systems involving queuing problems can be 
solved by simulation. In a queuing problem, an arrival 
occurs requiring a service to be performed. The system 
either performs the service or keeps the demand on a waiting 
list until it can execute it. PERT networks can also be 
Simulated. ALescheduleiofeevents on s"*critical path" hasito 
be found among a specific number of activities required in 


order to achieve a task. 


2.24 Applications in Ecology 

Becosystenst faseea Tawhole® areiihard Sto Bstudyteeveneny 
simulation because of the great number of variables to be 
included. Most of the simulations of ecological models in 
the 1960's involved small numbers of variables, and usually 
Gnby*,oneheor iitwo Species. stBut recently, because of social 
pressures and attention given to ecology and resources 
Management, ecologists have increasingly begun to simulate 


ecosystems. Examples of two approaches are given here. 


sacubt ata te! 30% snag siping: 4 retin aie —_ine Pe 
af apa eee ite sy alll ia adhd aware Aggespsasa yam 
Leworie | oF Nioiias0 Badecup. & oF wnoisntiiats yd bontom 
elie ye 97 babidtae Sada) soda (la paévayges satinoo a 
ugidigad’ » pd-tagieS ss sage Te bapties gig SaaGrive tsitde 
oi, waka» ago “ataodiec TATW | wat SSH IERS ase 24 Sevan Peak ay 
oe nal “ives lenititro" ip e&qavs) Sm plabedoe) é be to lneke an 
wi buries coitivitse' Se rdepn csphoeyzes sci Douek ed) © 
et £ overdose of aaa )/ 


yootow rk migksaabiqgh Sak 
74) seve) “Ente OF mul cette shade & He. -eneteyzand 8 of 
o¢ [e+ autdataay to 'zedpen tasty gat Xo, setec9n Ausra lpm 
Gh ptdQod Testeoloss Yo eaohrsteala sd To sen |.hetwiemh)” 
pitenye allah cic erie havieviir ease | 


19 


2.2.4.1 The St Kilda Ecosystem (Milner, 1972) 
SUPKiidatisa@analisland,; tpartebofa aAeinataonaletvature 
Reserve Min the vl Unitedtestates: It is a simple ecosystem 
since there are only a few vegetation types and only one 
major herbivore, the soay sheep. The model was built to 
study mainly the sheep population and the vegetation on 
Which they feed. The main steps of the simulation were: 
ia “Calculatetetie® growthwrof plants as’ a function’ of the 
temperature and the weight; 
2. Calculate the number of sheep as a function of body 
weight and age structure; 
3. Harvest from population of sheep; 
4, Calculate weight of plant material removed by sheep. 
These four steps simulated one day onthe island. The 
syst2m was intended to predict the future number of sheep, 
but did not do so very successfully. The conclusion was 
that either the model was structurally incorrect, or that 
one or more of the equations were wrong or had coefficients 
whose values were too far from reality. But still the model 
improved the understanding of the system and provided useful 


data. 


2.2.4.2 The Jabowa Model (Botkin, 1977) 

Trees compete for resources: light, water, minerals, 
and a place to grow. If a tree grows faster, it shades its 
neighbors who die without light. The growth, death, and 


reproduction of trees are functions of light, moisture, 


pres 20708 acm niwaie! 
oun yin rie” Say iuksaitaney: # ere 
od thsedh ul-thbde Sa speed “LAO _s sevontdand a, 
av anitsedapay ate fits doisetiqgog. grade a0 qfahn® i , 
‘Sey, ddtisl nate ei Do earte ‘chen age est pods tot je 

wip SH @erferwi wie. erdete 20) uiswedp .\at2 statyoles © 
<tip tow) oe. lune ee ee 

ghod) 20° -¢olrcrnt- ‘5 “Ga ojeete io edete sit .corloodgs Per 
| | [qapigia ge, spe bas |e a 7 

i 


tteeade io aottaloqog now teovIse- a 
aioete yd hevowas iczteten 361g Ve thGi aw! Ste Lewsss ah. 
ooo afinefead, edt ao” .y¥st ae fpeorvalretea “wots » tes vaae? | 
yeeede 930. asdulin Sty de* eit joLhety of Sebastian cae ayn 
saw poreoreaos atv cylietanepole Yaev oe oS Jen 646 ot 
jee 20) <dt@cdoun? ylieimsooo7. T Lshom giz asotie al 
senaioddsets bee To pyew wrew eqoirhips edd to esoe aad 
Poh ate ALize ao eyoruihsa wos! 225 ood aise captev | 


tulpan Aohivorq ice wareye ail? bh wiifantsisbap alt feve Q 


: oie vey) anos auodol a Tf 
te ene ao seanomnia’ x8? pone svat. 


time 


ee 


ae it om 


20 


temperature, and nutrients. Given equations for growth as 
functions of these variables, the Jabowa Model simulates the 
growth of individual trees on small forest plots. The 
program is interactive and uses a screen, a keyboard, and a 
Pughewepens The user can plant or log any number of trees 
from different species. Time in years can be specified 
along with the output interval. The program can also be run 
in "batch" mode for long runs. 

The model was tested and proven to work at least in the 
sense that it reproduced the general events. Simulation 
experiments were also undertaken to study the response of a 


forest to perturbations and manipulation. 
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CHAPTER 3 


SPATIAL DYNAMICS IN COMMUNITY ECOLOGY 


3.1 INTRODUCTION 

thesedensity?tof fa population#whashamportant effects on 
VEstegrowth: Overcrowded areas do not contain enough 
resources for their inhabitants and the death rate 
increases. In undercrowded regions, mating might not be 
sufficient to insure survival. Unless all individuals are 
distributed evenly, the density effect will vary from place 
to place. Spatial heterogeneity and ability of movement 
thus influence species interactions and dynamics. 

There are three patterns of dispersion in populations. 
Clumped or underdispersed populations occur when organisms 
are attracted to each other. The population is spaced or 
overdispersed if individuals repulse one another. Finally, 
the population can be randomly dispersed. 

Space Scanitibe itreated® un®laredisereter orticontinucius 
fashions S0rganishms can live in discrete habitable units, 
like the shoots of a tree. Or space can be considered as 
continuous where individuals occupy any point in space, like 
treescin’a forests) Finally; tis sometumes impossiblectto 
delimit individuals so that they can be counted. 

Studies of spatial patterns of populations of organisms 
have been popular in the last twenty years. Most of these 
studies, however, involve sessile or sedentary organisms and 


rarely has time be included as a variable. 
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Different approaches have been used to study spatial 
patterns in a single-species population. 

Before Lloyd (1967) defined "mean crowding", the usual 
measure of crowding was simply the average number of 
individuals. Actually, factors like the spatial pattern of 
distribution, or the way individuals react towards each 
other when they meet, or even the way they move, can 
influence*the degree of crowding. To include the «spatial 
aspectHrotne®theradsaspersioneofigaspopuhation; ebloydedckined ha 
measure of crowding as the mean number per individual of 
other individuals in the same quadrant. He also defined the 
measure of crowding when samples are used instead of 
complete data. 

Hatrston?, Hilieandarutten (1971), sanmplingt*ambieid ein 
Anh Arbor, Michigan, discovered "the failure of the sets of 
Ssampkeswitow reflect» thes distribution }coferthe ) linda vidual 
species ina satisfactory manner" (p. 343). But they also 
found that sampling in a grid pattern was far superior to 
Sampling at randon. 

For any motile species, patterns change with time, 
making "dispersal" or "diffusion" studies more appropriate. 
Dispersalt is "the apparently aimless, underdirected 
movements of animals that seem to be wholly random" (Pielou, 
1969). Organisms disperse to find a suitable habitat. A 
population responds to local variation in the environment 
with movement of individuals between areas. 


Holgate (1971), in his survey of results in random walk 
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theory, elaborated three models according to whether 
IGdLVidualsihave noskattracti#ony tokuafahomes@baseresoperthe 
attraction diminishes or increases with distance. Levin 
(1976) provided a general framework for modeling the 
dynamics of spatially distributed populations. He used a 
continuous time approach but the environment was subdivided 
into discrete islands or patches. The species density was 
then defined in terms of differential equations, including 
the intrinsic rate of increase of the species, the carrying 


capacity of the environment and a coefficient of diffusion. 


3.2 MULTI-SPECIES SPATIAL INTERACTIONS. 

If many species live in the same environment, or if 
they have some effect on one another, their patterns may not 
be independent (Pielou, 1969). By statistical analysis of a 
table whose entries specify the number of cells or units 
that contain individuals of every species, one can determine 
if two species are distributed independently of each other. 

Murdoch (1977) studied the stabilizing effects of 
“Switching", where a predator divides its time between two 
different prey areas, spending more searching time in the 
area that contains more prey or where the capture is more 
"rewarding". He came up with a model which describes the 
number of a given species eaten per unit time as a function 
of the transit time between the two areas, the time taken to 
"handle" the average individual of this species, the density 


of the prey species, the successful searching rate for the 
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prey species and the search time. Results fron experiments 
with a predatory fish (guppies), a surface prey (Drosophila) 
andfasbottomiepreye (tubificid 4vorms) everifiedt thei) basic 
assumptions of the switching model. 

Investigations using simulation have also been done to 
Study spatial behavior in predator-prey systems. The work 
of Huffaker provided the experimental frame and data for 


many of the studies. 


3.2.1 Huffaker (1956,1958) 

Huffaker studied the use of predators as a density 
control tool. Strawberry crops in California suffered heavy 
damage from the cyclamen mite Tarsonemus pallidus Banks 


unless a predatory mite, Typhlodromus cucumeris Qudemans, 
waS present. Both field and greenhouse experiments were 
conducted. 

A female cyclamen mite lays about three eggs per day 
foun EoguE tons Livestedays: Only about five per cent of the 
aduitapopuldtion sisemales/ eByecontrolbings ,its ire productive 
rate and its searching behavior, the predator can both hold 
the densityhofivethe eprey tat Slevels@c<nofbnjurioussttoe tite 
strawberry plant and also can survive at low prey density. 
The predator usually delays its appearance anvailfieladt until 
the prey is well established. 

The main goal of the field and greenhouse experiments 
was to find a way of controlling the predator population. 


Sampling was used to estimate the different populations. 
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The conclusion was that the predators "exert very effective 
control of the cyclamen mite in third- and fourth-year 
fields in this area when their activities are not inhibited 
by detrimental chemical treatments used for control of other 
pests"! (Huffaker sands Kennett, 7:1956,, )p. 198): Forteathe 
greenhouse tests, seventy-two potted plants were arranged in 
two randomized groups. One group was kept predator-free and 
the corresponding plants*rapidly d#edaplGhetresnults fromthe 
studies suggested that deliberate infestation of new 
plantings by the predatory mite could be a means to achieve 
equilibrium and preserve the crops. 

The objective of the 1958 experiments was to create an 
ecosystem where the predator and the prey would continue 
living together so that they could be studied in detail. 
The prey was the six-spotted mite, Botetranychus 
sexmaculatus and the predator, Typhlodromus occidentalis 


Fach universe consisted of a metal tray containing a ten by 
four array of oranges and/or rubber balls. The predators 
and prey were free to move onto or leave oranges or balls 
but not to leave or enter the universe. The utilization of 
rubber balls made it possible to change the total food 
substrate (oranges) and the degree of dispersion of that 
substrate. At intervals of eleven days, one quarter of the 
oranges (the oldest) were replaced. Sampling was used _ to 
count populations. Results were obtained both when the 


plant food was readily accessible, all in one location, and 
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when the food was widely dispersed with intervening rubber 


balls. Groups of prey with and without predators were 


included. Different areas of food were established by 
covering part of the oranges with paraffin and paper. Some 
of the results are summarized in Table 3.1. "Adjacent food" 


means®’a S2-orange, feeding *areamion» 14) adjacent“oranges. 
"Dispersed" implies the same feeding area but dispersed in 
the tray. The last type of experiment used a 2-orange 
feeding area on 20 oranges dispersed in the tray. 
Conclusions from these studies were that the 
environmental area must be large enough to permit co- 
existence of the predator and the prey, yet not so large as 
to hinit interchange of the biotied participants: 
"Migration" (or, here, simple movement) is necessary to 
avoid extinction of either the prey or the predator. It is 
possible to establish complex universes, utilizing wide 
dispersion.and also including partial barriers, in whichrthe 
predator-prey coaction would not be self-annihilative. 


Table 3.1. Mean prey population per orange 
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Predator 
Food absent present 
adjacent 4700 250 
dispersed 3500 4056 


alternating with balls 3300 1800 


a ek cas) cous Ca es ca es wees a ee aan aS tse a a Se ae eS SD GED AO SOTO Oe EP ES OE Sm ae, Se A OS Ee ee ee 


Phogt seenaeae® oat tee! sti: ne Sy ee | 


ees Dienst te ao wees ERGest agtarp~S 6 eéuom, jl 


<i pebeaah Mt Bere eeenec? dase S@y @abighe “Soacsqgeee 
eunadiant w bein Stenigegrye te Sats (feet ii? .9a29, ee 
‘vets dar ai So@nedeal) apmprco0t ao eee Galion 

eat a 37° aapbuve spent so73 anoszeione).— ‘4 
ce Yiweed oF dptone Soael se ?20@ sazm Levanesonteae. 
a ib! a6 toa Soy 5192q) 502 Dap 1otpbegh, was to eon sieite 

sae eanisieitg.-xs totic cal? to siadp wy res jinti a2 
ad “teeeahar 3k “faneetes Gigdie 0380. (3e)° “eotseamel? 
at at “sefabous 944%. 90 Yady alyoaes tis ry Aelruans» Liovs 
olive peisélljye veteievkinn teiqnis aeiLagte> of efpaneng 
J4 eplde bi vee Perot deotier poisilow! Oels ‘fre soleregEem, 
seveteal f¢hace-3 fee ed for Pipew nossneso 1929-20 10eeTe ; 


“a gins tials 


7? 
1a2hh 17 af oe 
c... mc a 
-_ MARES ePNRE AN RRS ARISES MER YEE BRE 
al =~ ee ee eee o= 5 oe <= ae -——<2- == ‘ae < 


is halal act ea a _ 


Pa | 


3.2.2 Murdie (1971) 

Murdie suggested that the spatial distribution of prey 
dnfluencesethesability of!predators to “search "successfully 
for food, and that adaptive behavior is shown by the 
attacker. A random walk search by the predator was tested 
with different prey density and spatial distribution using a 
computer-Simulated model. 

Some of the models describing predator-prey 
relationships include a coefficient of "killing efficiency", 
Which might reflect the competition between predators or the 
searching ability of the attacker. For example, some 
predators change their pattern of searching when they 
encounter a prey. An intensive search around the discovery 
area begins#! whichescan abeelirurttulseifdwethesgpcrey stare 
aggregated but detrimental when they are dispersed. 

The environment is represented by 150 quadrats of 10x10 
cells each. A guadrat is a square area of vegetation marked 


off for study. Using the negative binomial distribution 


= . xXx re = )! 
7 =ea Geena) eat 


where X¥ is the mean density per quadrat, a frequency is 
determined (f, = Npy, where N 1s the number of quadrats). 
Different k and x values are used and a number of quadrats 
| £ | are each allocated x prey. The prey are either 
distributed at random within a quadrat or a random location 
is chosen in the quadrat around which the prey are clumped. 


Only one prey is allowed per cell. A fixed number of 200 
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predators are used, each of them being assigned to one 
quadrat at random. All predator-prey encounters resulted in 
a predation event. After twenty successive failures, the 
pEedatorsis) assumed “totadie. Foret soney of Pt thet oruns eta 
satiation level is set at ten prey. 

Results showed that the clumped model has a lower level 
of predation then the non clumped model, thus suggesting 
that aggregation has some protective value against random 
walk predators. The introduction of a satiation level 
Significantly modified results. As expected, the number of 
prey killed waS proportional to the number of prey 
available. In the record of the discussions following the 
CONUCELEDULILONP ORNS ES Waters from the U.S. Forest Service 
notes that prey movement should be included and Murdie 
himself suggests the introduction of reproduction for both 


predators and prey. 


3.2.3 Maynard Smith (1974) 

Migration is believed to have a stabilizing effect on a 
predator-prey system. Maynard Smith developed a model ‘based 
on the following assumptions: 

inivhethabitat“was* divided sintonaticells"pawhich) couldsobe 
considered either as a continuum or as a discontinuous 
habitat; 

2. Then effects soot immigration and emigration were 
immediate; 


3. Migration was "conservative", in that the system was 
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closed; 
4. The environment was uniform. 
By graphical analysis, Maynard Smith showed that migration 
cannot stabilize an unstable predator-prey system to insure 
coexistence. A model was built where a cell was considered 
as @an Stlonange" (see,Huffaker7ei958) sfetThescell could be in 
eight different states: empty (E), containing few prey (HA), 
increasing prey (HB), many prey (HC), many prey with few 
predators (MA), many prey with increasing predators (MB), 
many prey with many predators (MC), and few prey with many 
predators (MD). Probabilities were defined for state 
transitions. The general conclusion, after the behavior of 
the model was analyzed by computer simulation, was that 
persistent coexistence of predator and prey could easily be 
achieved. Conditions which helped this coexistence included 
high level of migration and presence of refuge for prey, 
restriction of the migration period for the predators and a 


large number of cells. 


3.2.4 Ziegler (1977) 

Most analysesWiot ‘spatzal ceffectsy Pav tecology have 
employed the assumption of continuous population exchange 
mechanisms. Ziegler considers a system where migration can 
occur only at specific points in the predation cycle, such 
as when the prey or predator population of a cell reaches 
the carrying capacity of that cell (like the Maynard-Smith 


model). The results of his simulation are given here. 
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The base model (see section 2.2.3) is specified by the 


Lotka-Volterra equations 


ax7dat 


Ul 


AX#(1Re ex): SOx y 


dy/dtt= -byit xy 


where 
x = number of prey, y = number of predators 
a = growth rate of the prey 


b = death rate of the predator. 
Paraneterst*i| and 602, 9®the>ipreystand tapred ator teextinction 
levels are defined, which encode the notion of "sufficiently 
large" number of organisms. A cell (lattice point of an 
infinite two-dimensional grid) is then defined to be locally 
Vunstabiemeretba<thitandtfal<iL27(d=b) > 
The assumptions of the lumped model are: 

1, Each cell goes through the following sequence of states: 
EMPTY, PREY (containing prey), PREY' (containing maxinum 
number of prey), PRED (containing both prey and 
predators), PRED! (containing maximum Humber of 
predators), EMPTY. 

Paaiiaratiouwot tprey #fron feel IMistokcel te jrocctrstoniy Wien 
eal] weiss binePREY Sastatesand@celing asteMPry: 

BseMigration of predators: fromycell 1 toscell®} occursvonly 
when cell i is in PRED*' state and cell j is either in 
PREYSocrePREYUastate:, 

4. Migration from cell i to cell j has no reverse influence 
on cell i. 


When a cell enters the state PREY', the neighboring cells 
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are scanned and each EMPTY cell is converted into a PREY 
celtepwithetprobabitityseths In the same way, when a cell 
enters *CheestatesPREDY) the adyjacenticelissin BREYtor PREY! 
State are converted to PRED state with probability p. 
Simulations were run using different values for p and h and 
Lesultseiwveres gidentical:'to Maynard=Smithise (sees section 
3.2.3). The mathematical model accounting for the results 
took the form of the Lotka-Volterra system and could explain 
the ability of an unstable system to stabilize through a 
discrete population exchange mechanism. 

An ecosystem is usually complex and a single model to 
represent it, however desirable, is not practical. Ziegler 
(1978) suggested that partial models be built corresponding 
to different experimental frames and that the computer be 
used to integrate this collection of models into a 
synergistic whole. These models can be partially ordered by 
use of morphism relations which induce "a mapping from the 
parameter assignments of the finer model to those of the 
coarse ones" (Ziegler, 1978). AS an example, Zeigler 
considered the universe employed by Huffaker (1958). Four 
related models were built representing different refinement 
levels. The first (base) model was of stochastic 
differential equation type and represented one local state 
(situation on each patch). The second (lumped) model, of 
stochastic discrete event type, kept track of the same 
variables as the first model but updated them only at "event 


times", The third model ("occupancy") is an abstraction of 
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the second model in which a small number of states was used 
to represent the local situation. The fourth model ("random 
phase-space") described the global behavior of the third 
model. 

All the models could qualitatively explain the observed 
persistence in Huffaker's universes and the "more refined 
models were able to give close quantitative agreement as 


well" (Ziegler, 1978). 


CHAPTER 4 


ECOSIM 


4.1 INTRODUCTION 


A system called ECOSIM has been developed to 
interactively simulate ecological phenomena on a digital 
computer. ECOSIM allows a user to model and simulate the 
behavior of a community of many species in both space and 
time. The program is capable of handling large "ecosystems" 
including numerous species and many individuals. Easy to 
operate, ECOSIM's usage requires no knowledge of computers. 
A dialogue is established with the user to allow him to 
build his model. For each species, questions are asked 
about its different characteristics and parameter values. 
Choices are given for output intervals. Simulations can be 
run on previously saved or newly created models. Editing 
facilities are provided so that between runs, the user can 
alter the environment in any desired way. The main features 
and options of the system are outined below. 

Hexagonal tesselation. The environment is represented 
as an unbounded two-dimensional grid of hexagonal cells. 
The hexagonal element is preferable to a square or triangle 
because it has six unambiguous neighbors and permits better 
approximation to arbitrary delineations. Rach s cain can 


contained? Teeortinoneatindividualsieand/or 0) “loeor, nore 


barriers (to be defined later). 


33 


os ‘uqeniies soot ost onngo7d | Bites date Poo a 
jetieih. «po. anvapmedt (eatpolese etelowke ylevksossutal 
atv afeiare 24) Leow of sseu 8 | evolin MERQUE  .aechqued” 
hue soace vod ov ssiteqp yres 30 Yieauaoo a 20 toiveded 
thindeurmone” soret taifinad to eld@eqe =e Bempesqg sd? owed. 
St - Weak al swtivihol vieaour ealinsge éenmsena pa thuiaat ; 
asseqeos 46! GW 1vénd op ea Dipset Bye) eheseos ( etrmega 
of wid wetde! oo Yeen afd) ARTY haddildaves ef snpoberd Lon 
hedue ‘sxe gvoiseaip “satanye dose 20% sfeboe aid Miod 
enn Lev tse enn ae “bra entip ied osiedo sapuetaeb es seede 
at agwoudstisvelupi? salivrernd sone 307 osvtp 238 eestos? 
pxtarép .efebon haseers Wie 37. Omvas Yieieivetq # —_ 
as5 eeu ml? .anns ohevdei je). oe Bablgoiy #26 setnationy 
saiutes? nike eft .yov fealweb yrs Wh aismnoives of7 299 ka a 
Wolsd Roal wo es wagaya oft 20 aaniege are 
hetacséryov at taseodives ody -Heddkleetus Lamuounel. ~ 
efpaciat 20 spaupe a od slasictesy 22 trenets |b propped sad? 
geste -wtieing fine Bz0dtpion auemp romano vis cee i ounesed 
asa diea sine ere eee ox sod | 
‘aie wt ve a 
7 an ry Winiten.9¢ 0% wis 
‘ 
if 


Sian a | 


34 


Species and individuals. Any number of species can be 
included in the model. Any number of individuals can belong 
foGeachespecies. The only practical limitation on the size 
of the simulated ecosystem is imposed by the availability 
and cost of computer resources. 

Died WSt orp ton. After specifying the initial 

number of individuals for a species, the user either enters 
their locations (as coordinates) or selects parameters for a 
stochastic distribution (uniform, normal, exponential, or 
Poisson). 
Processing level. A processing level is associated 
with each species. Processing levels determine the order in 
which species are processed during a time step, individuals 
of the motile species with the lowest processing level 
moving first. 

Mobility. A mobility factor determines the probability 
of each individual of a (motile) species of moving in any 
time step. The mobility factor can be dependent on the 
hunger level if the user wishes. If so, parameters of the 
quadratic equation for the curve passing through (MIN,1), 
(MED,MOB), and (MAX,0) are determined, where MIN is the 
minimum satiety level (to be defined later), MAX is the 
maximum satiety level, MED is the middle point between then, 
andusnoOB@sis ithe’ imobvirty Goetiticirent. The mobility 
Coefficients corresponding “too any satiety > slevel scan = be 
calculated from the guadratic equation. 
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factor determines the probabslity phofy.e li ninatiugeean 
Syavviduat@etrying@ tor occupy isthetesane® celiatas Panother 
individualeeof i thei tsame!l ispeciesoerifythéatser Gspecitieska 
random walk, the occupied cell and each of its neighbors are 
assigned a probability of 1/7. Alternatively, the user may 
specify the probability of movement in each direction. 

Predator-prey relationship. For each prey species, the 
user assigns to each predator species a probability that a 
predator will kill a prey, as well as the maximum number of 
prey ajphedator) canjkiblepertunitetime. 

Breeding specifications. For any species which is a 
predator and produces offspring, the user sets a breeding 
interval (in time units), a probability of breeding, and the 
mean and variance of the number of normally distributed 
offspring. At the expiration of the breeding interval, some 
new individuals may be "created" in accord with the above 
parameters. The offspring are randomly distributed in cells 
surrounding the parent. Each newborn is given the same 
satiety level (see below) as the parent, a breeding clock 
tinertiote zéro,e and@ thesaspecies' trophicetlevel ponobnlity, 
predator, breeding, satiety and exclusion characteristics. 

Satiety specifications.) The’satietysofesan individual 
is indicated by an integer value ranging froma minimun 
value indicating the starvation point (the individual dies) 
to a topt ivalueisindicating Sthe)"satiationy levelia(tie 


individual does not eat even if he encounters prey). An 


increment value chosen by the user is added to the satiety 
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level each time the predator kills a prey. A decrement 
value is subtracted for each time step the predator does not 
f£indleayspreyotinesits cellietForsveach predatorespecies;ithe 
user sets the mean and variance for the initial (normal) 
distribution of satiety values. 

Exclusions. The user may specify “incompatible" 
species. "If species*i éxcludes species 4, an individual’ of 
species j will not be allowed to move into a cell occupied 
by an individual of species i. 

Global modifications. At the end of each time step, a 
general function (specified by the user) can be applied to 
an existing species population to change the number of 
individuals. The global modification is specified as a 
Fortran expression of the form f(x,t) where x is the actual 
population, and t the actual time step. For example, if the 
user specifies 2*x*sin(t) as a global function for a given 
species, the population of this species will double at every 
time step, weighted by a factor of sin(t). EH TX (et 1) oa 
x(t); deletions are made at random. BECEXEed) Sax CEG 
additions of new individuals are made according to the 
initial distribution. TEATthe® tnatialeedistuibutiond was 
ex pressed aseespecific) Rcoordinates ms etheauuser will be 
requested to enter as many coordinates as required for the 
"ew"! individuals. A satiety level corresponding to the 
satiety initial distribution is given to each individual. 


This feature can be used to simulate migration and seasonal 


variation. 
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Barriers. Barriers can be established either to create 
refuges for. certain species’ or to impede or prevent 
mobility. Fors eachéebarrienasgtype, art he user provides 
peobabilt@ties tof® crossingetiorudany number ofespeciese(the 
detautteevalue’ Ofor®ithewaprohabilityy ofpecrossinge@eist 13.2 
Barriers can be distributed either at random (uniforn, 
normal, exponential or Poisson), by specifying coordinates, 
or by giving the four side coordinates of a rectangle. 

Output capabilities. The user can have the environment 
displayed graphically attihkanyreterminals Displays are 
provided initially, at specified time steps, and at the end 
of the run. The environment is divided into windows of up 
to 16x16 cell regions, and the user can display all regions 
or only desired regions. 

Editing. The user can alter the ecosystem 
specifications before initiating or continuing simulation. 
Species can be created, individuals added to or deleted from 
existing species, all species characteristics modified, and 
individuals moved from their location to new cells. In the 
same way, new barriers can be created, existing barriers 


moved or deleted, and new barrier specifications defined. 
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Savi 
actual simulation, the existing state of the model can be 
saved in a file whose name is supplied by the user. in ra 
subsequent run, the user can request that the simulation be 


run using the information on the environment contained in 
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Run. The user supplies the number of time steps, the 
output interval, and any stopping conditions (if a given 
species reaches a certain maximum or minimum value). The 
Simulation is then carried out, using either a newly created 
environment or one that was previously saved, and output is 
produced as required. 

Dialogue. The data structure needed for the simulation 
is built while the user is answering questions from the 
program. Error recovery is possible from almost any part of 
the system. The answers of the user can either be spelled 
out {like YES .or..BARRIER).-or.abbreviated, with the ’first 
letter (Yeon rs): When an invalid answer is given, the 
question is usually repeated. Sometimes an explanation is 
guyven Oftthe. error,-ifs«.the, Situation. .iSs.n0t,.clear: At 
certain points, a null answer has to be given to terminate a 
cycle of questions. Asetulleadescription,otfythe riles is 
given in the User's Manual (see Appendix). A hypothetical 
dialogue f£oldows -ieUTlustratinge.vartouss.optaons WoL ~ the 


system. 
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4.2 System Structures 


The system is composed of eight main modules including 


the MONITOR. A schematic diagram is shown in Fig. 4.1. 
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Fig. 4.1 Schematic representation of ECOSIN. 
The set of records representing the data for species, 
characteristics is 


barriers, predation, and exclusion 


accessible from most of the modules. Bach \enodulewtis 


considered below. 
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4.2.1 MONITOR 


The MONITOR is the control center of the simulation 
system. Initialisation and contreotenofteabiteginternal 
Operations are done in the MONITOR. The main sequence of 
operations is: 

1. Explain 
Ze Create a model or get it from a file 
Sagaoutpar 
ioe Edit; save, *quit; or run 
Sse Legedit edit.andigoeta ks. 
6. If save, save and go to 4. 
Le BPLLeVrun;y Luneandeqopcons: 
So PERE GuUaty cure 
QSOIEdiG, Saveyuq ult; cor continue 
10. 88LG; e€drt,tedit andegottons. 
THiS aelicsave,~savegand gotta 9: 
2S melt run, fruneandggqosgto ha. 
P3RPLELeCuUst peg uate 
AGingdel caneithustrbe rsaved  yhbeforeveand#aiterythepactual 
Simulation, allowing study of the model before further 
Simulation is attempted. All steps (except 2) are optional. 
Once the user has examined the results of a simulation, he 
can simply restart the system, obtain the model froma _ file 
and use it to "continue" the simulation. The random number 
generator is initialized by the MONITOR using the computer 
clock as a seed. Also initialized are the default minimum x 


and y value and maximum x and y value. Record definitions 
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aresgiven. for species,.individuals, barriers, barrier cells, 
predation, exclusion, species/barrier and halting 
characteristics. Pointers to the different records are set 
to null. Tables are created for the parsing of the global 


functions. 


4.2.2 INPUT/OUTPUT 

Several procedures are included as input/output 
facilities. EXPLAIN is used to refresh the user's memory of 
the facilities of ECOSIM. Since the information provided by 
EXPLAIN will not be sufficient for a new user, a user's 
manual is provided (see the Appendix). 

The procedure CREATE initializes the main parameters 
and controls the "conversation" about the model between the 
system and the user. Questions are asked about the 
characteristics of the model being built and data structures 
are created. The procedure can be called from the MONITOR, 
or from EDIT if a new species or barrier is introduced in an 
existing model. Error recovery is done simply by repeating 
the guestion if the answer is invalid (like a real number 
for an integer or an illegal ‘character in) “a number)- 
Breeding and hunger characteristics are requested only if 
the species is a predator. A normal distribution N(y,c), 
where Mis the offspring mean ands is the offspring 
variance is used to compute the number of offspring. If the 
mobility is dependent on the satiety, the coefficients of 


the quadratic equation corresponding to the given variables 
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(see Fig. 4.2) are calculated and stored. Ata given time 
step, the individual is then moved if a number picked at 
random is less than aHUN2 + bHUN + c where HUN is the 


satiety level of sthe individualvat this time step. 


Wake aHUNBOT2 + bHUNBOT + c = 1 
| aHUNTOP2 + bHUNTOP + c = 0 
i\ aHUNMED2 + bHUNMED + c = MOBCOEFF. 
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To generate the individuals, the user first enters the 
number to be created and if they are to be placed at random, 
me not, the given number of individual locations is 
requested. The user can either enter the locations one by 
one or retrieve them from a file. If the individuals are to 
be placed at random, a distribution and parameters are 
chosen by the user. For each individual, two uniformally 
distributed random numbers (one for each coordinate) are 


transformed by standard procedure into numbers with the 


given distribution (Fishman, 1973). 
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If the user wishes to use a model already saved in a 
file, the procedure GETFILE is called which handles the 
creation of data structures and transfers the data from the 
file to the records. A check is*made to see if the file 
exists and if not, the user is prompted for another 
filename. 

Before or after each simulation run, the existing model 
can be stored ina file via SAVEFILE. The procedure goes 
through the data records and saves all needed information 
for subsequent runs. Values of important variables values 
like the seed for the random number generator, the number of 
species and barriers and the minimum x and y and maximum x 
and y are also stored. The hashing table is not stored 
Since hashing is done on newly created records or on records 
retrieved froma file. Other procedures are used to output 
information about the different species and barriers, if 
requested. 

At given output intervals, the present state of the 
model is displayed as windows of up to 16x16 cells in a 
"grid" form. Individuals and barriers are shown with their 
display symbol. At the beginning and at the end of each 
run, and also at output intervals if requested, statistics 
are given about the number of individuals left of every 
species and their characteristics. 

Throughout the simulation run, results are saved in a 
temporary file created at the beginning of the run. Before 


guiting, the user can request that this file be made 
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permanent and/or printed on a batch printer. If the file is 
to be made permanent, a filename is requested and the file 
is stored permanently if enough space is available. 

A plot option allows the user to plot the population of 
desired ispecies atiweeach o.time) ostep: For each species, 
population versus time step is plotted on the Calcomp 


plotter. Each answer YES to PLOT? generates a new graph. 


8.2.3 EDIT 

Most of the model's parameters can be altered via EDIT. 
New Species or barriers can be created or existing ones 
modified or deleted. Any general parameter of a species or 
a barrier can be altered. A species, for example, can 
become a predator or can cease to be one. Individuals or 
barrier cells can be moved or their characteristics changed. 
Predation, exclusion, and species/barrier parameters can 


also be modified. 


4,2.4 HASH 

Frequently during simulation, the contents of a given 
cell have to be accessed to determine what individuals or 
barriers it contains. A hash table is built so that, given 
the coordinates of a cell, a pointer runs through the data 
structure to every individual or barrier cell occupying the 
given cell (see Fig. 4.6). Before hashing is done, the size 
of the table is determined: the square of the longest side 


or 100 (whichever is less) is assigned as beheigsizew@ of the 
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table. For example, if the environment is a 5x3 cell grid, 
Lucrstzesororhe: hash Ptable twit iebe §25¢ 162). Then the 
hashing function assigns an integer number to the cell 
(modulo the*sizej of? thel’table) according tolatsS position in 
the environment (see Fig. 4.3). For a given coordinate 
pair, the corners of the "square" to which it belongs are 
found. For example (from Fig. 4.3) the corners for the cell 
Giant) (whach?acorresponds =to, cell916) \areut0,4¢947549.8, and 
22 {second square from the center). To find the corners, 
the maximum distance DIST of the cell from the center cell, 
i.e. the cell in the middle of the grid, is calculated. The 


four corners are then 


CORNER1 


(DSH) *4 (DASia st 1) eee jee 
CORNER2 = CORNER1 + DIST * 2 
CORNERS#= CORNERZ! 4 DISTAs 2 


CORNERS# +2 Di STSt 2 


CORNER4Y 
binelly;,.according to" athe § position (ob the cell, in the 
Square, its distance from the appropriate corner is found 
which then added to the corner itself gives its hashing key. 
If the number is greater than the size of the table, it is 
divided by the size of the table and the modulo is taken. 

A pointer goes from the hash table element to the first 
individual of the corresponding cell, and then from this 
individual to a next individual or barrier also occupying 
this =cellt ‘or |‘ haShing “tor *the »same. table® position. «The 
hashing function is used each time an attempt is made to 


move an individual, to see if there is a barrier the 
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Fige 4.3 Example of the hashing function assignment. 


Pudividualgecannot s#eross) ¥oresifidgexciusiont ongpecompetition 
prevents the individual from moving. After hashing is 
performed, a check is made to (verify if any individual 
occupies a cell containing a barrier that he is not supposed 
toeceoss. (With probability 1.2.0). Lf ‘so,- ay warning Fis 
issued, indicating the conflict and the user is returned to 


EDIT to make the appropriate changes if desired. 


4.2.5 COMPUTE 

For each species, the user can specify a function of 
Fhewecureent!) Sspeciesiasizestandiecurrentettine steppetothe 
applied at the end of every time _ step. RAGrortran® Like 
expression has'to be entered, as a function of x and t, like 
K*¥*2 or x*(t+1). If the result is negative, deletions are 
made at random: If the population is’ to increase; additions 
of new individuals are made according to the initial 
distribution. COMPUTE is used first from CREATE to analyze 


the Fortran expression for correctness. After each time 


otremdezuss golsiuwr poddend ghz — Ene pa 


toltieajevs 20 solvirdowe |7f ae apome _jotess Isehtedbatl - - 
#l. peidest Garis .polivon wort Levdkvitas edt aragiaig — 
bawlgy fhek Gok 12 BPicev (od Shee ek apaia @  sheutotIeg 


igtneyye cog 22° oR) TAR - *weteRe a: vodeioveso bien 6 eoigeaeo” | 7 
af oniece® © G2. TF eto. VrsALSaUIN, its#) sa073 oF 


i 
+9. Sapyetet a2 Tee! of Bob sr FT Qos, ost pakraw last sosaeat : 


euatech 42 @ewieds ‘eratsqozqge, ahs shee o3 PEGS, 

. —- 

FRUAOD eRe. 

ty «pobenwir? + ('¥oltegeeaeo ana sit. wellpage dase 20% ¢ Of -_ 

of e2 .gete eel? “tredat Bab inate) aetasge) faapree ae 
oii anette’ A quate swt yaave fe fos ade to Besiqae 

etki \¢ lon * Yo Wotan? n se ,heteene eh Od, and eoteseaqxe 

ete adoitaleh javlappen si: rivera sth MT «(Peeps 10 Ree 

esorrigse osenortnd of ei noitgatuqog sae 2T Gwebnee ae 

Katrids ots oF pipes. shew sib alsubivibas. van 


sittin on 37409 bad pense = 


53 


step COMPUTE is called to execute the Fortran expression and 


determine how many individuals should be added or deleted. 


4.2.6 RUN 
RUN is the main simulation module, controlling the 

sequence of events in the -systen. AIGiist Us first 
established to define the order in which the species will be 
moved (by processing level). The number of time steps, 
output interval and halting conditions are then requested. 
For each species, an individual is found, at random, in the 
list (see Fig. 4.4) and a second random number determines 
the a“direction" in twhicht etheslistsoitindividuals will be 
scanned. For each individual, the following sequence is 
executed: 

1. Adjust the age of the individual 

2. Move the individual if required 

36 Let the individual cat, if required 

4, Let the individual breed, if required. 
After all individuals of a given species are processed, the 
guoba la tanetionsofathisespecies isy applicdgrri situpexzistse 
Aftereall species»are. processed,a.check isi made \toeseelifia 
species is out of bounds, 1a. 82f vthesnumber of Sindtviduais 
Opea Spectestis@ibelowt themeinieniorsalinitehore®abover the 
superior jelimito tset | l\byethetusers {eif so; the ysinulationwis 
stopped and control is returned to the MONITOR. The output 


interval is then checked, and output is produced if 


required. 
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4.2.6.1 MOVE 

The cmobility coefficivent of the individual is 
determined according to the mobility variables a, b, c (see 
4.2.2), the hunger level, and whether or not nob Tit ye ais ita 
function of hunger. A random number is picked up. If the 
individual is allowed to move, the following steps are 
executed: 

1. Find in which cell the individual will move, according 
to the probabilities of the neighboring cells; 

2. Check this cell to see if the individual is excluded, or 
if competition or a barrier prevents him from moving in 
thismaice i): He then dies and procedure KILL is used to 
delete the given individual from the data record; 

3. If nothing impedes movement, move the individual by 
adjusting its coordinates and adjust minimum and maximum 


Xtandwyetiferequired. 


4.2.6.2 EAT 

Limthe® predator *'isythungry “(satiety Wevel < satiety 
top), a check is made to see if any of the prey occupies the 
sane Gell) as ithe predator. Ef not; “the 'sdtiety level ofthe 
predator is decremented. If a prey occupies the cell and 
tie sepredatommdd ids not" reach at lem tsataation (evel gor the 
maximum number:‘of kills per unit time, the prey is eaten and 
the check in the cell continues until the predator has 


reached the satiation level or has killed the maximum number 


of prey. 
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4.2.6.3 BREED 

According to the breeding probability, the number of 
offspring is first determined by generating a random number 
from a distribution N(u,s) where mM is the offspring mean and 
o the offspring variance. For each of the offspring, a cell 
is picked at random around the parent and if competition, 
barrier, or exclusion characteristics of the cells occupants 


allowsiit;athe offspring is “createdfinithis cebis 


4.3 DATA STRUCTURES 

The data structures are intended to save space and 
allow easy access and have been designed in accordance with 
the following constraints: species can move and grow without 
bound; editing can insert or delete whole components of a 
model in the middle of a simulation run; all members of a 
species must be found and processed at the same time; all 
organisms and objects in a given cell must be identified to 
determine if an organism can move there. Since only the 
last constraint suggests an array type organization, linked 
allocation of storage is the method of choice. 

When a new species is introduced by the user, a_ record 
(called the SPECIES HEAD NODE) is created containing the 
general information pertinent to that species (see Fig. 4.4 
andes. /imeen pomnter@allowsh docess® toeithenelrrst (especies 
record anda pointer from species to species permits access 
to*®eachedéf them-ashrequired: \"individualsysof each species 
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individual node (IND. CELL NODE) carries the organism's 
location (coordinates), age, satiety level, and number of 
killed prey. A pointer to the species head node allows easy 
access to the general information about the species. The 
doubly circular linking allows control of systematic effects 
which might arise from always processing the individuals in 
the same order. Each time step, a random individual is 
chosen and the list is traversed in a randomly chosen 
direction, so that on the average each individual is 
processed before and after each other individual equally 
often. 

Barriers are represented as singly linked location 
nodes (BARRIER CELL NODE) with pointers to a head node 
(BARRIER HEAD NODE) containing parametric barrier 
characteristics (see Fig. 4.5). A top node in the system 
contains pointers to the first nodes in the list which links 
the species head nodes and the barrier head nodes. 
Interspecies exc lusion factors (EXCLUSION) , predator 
characteristics (PREDATION), species/barrier specifications 
(SPEC/BARR) and stopping conditions (STOP COND) are 
maintained in separate lists. 

The last data structure element is the hash table whose 
nodes head the linked lists of all organism and barrier 
nodes whose locations hash to that table entry. (see Fig. 
i 6) Thus the last constraint mentioned above is 
satisfied. The hashing function is used “each time an 


attempt is made to move an individual, to see if there is a 
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CHAPTER 5 


SYSTEM EVALUATION 


Se 1 APPLICATION: MURDIE'S MODEL 


5.1.1 Background 

ECOSIM is flexible enough to allow study of interesting 
ecosystems. Simulation runs can be carried out step by 
step, introducing features of the ecosystem gradually, to 
verify their effect on the model. The system studied here 
involves Murdie's (1971) experiments. He suggested that the 
Spatial distribution of food or prey influences the ability 
of parasites and predators to successfully search for 
victims, and that predators show adaptive behavior in their 
searching patterns (see Section 3.2.2). After simulating an 
environment containing both clumped and non-clumped prey, 
Murdie concluded that clumped prey had a lower rate of 
predation than prey distributed at random, when predators 
moved by random walk. 

Murdie distributed prey and predators in an environment 
represented by 150 independent quadrats of 10x10 cells each. 
All prey encountered were eaten except in certain runs where 
a satiation level of 10 prey was set. After twenty 
successive time steps without encountering a prey, a 


predator would die of starvation. 
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5.1.2 The Model 

The model environment built using FCOSIM consisted of a 
10x19 cell region corresponding to one of Murdie's quadrats. 
Various prey and predator densities were used. The prey 
were either distributed uniformily or clumped using a Normal 
distribution. Most of the time, the clump was placed in the 
middle of the grid, so that the "boxing-in" effect would be 


minimum at the physical limits of the region. The predators 


were distributed uniformily. A barrier rectangle was 
established aroundev thes, regrvonetsor ithatiebothneprey sand 
predators remained within the 10x10 area. To Simulate 


Murdie's experiments, the prey were kept immobile while the 
predators moved by random walk (see Table 5.1, model type 
(T)ipe For some runs, a Satiation level was set at 20 prey 
(Table 5.1, model type (2)), with starvation after 20 time 
steps without encountering a prey. Other features were then 
introduced like mobility of prey (type (3)), breeding 
characteristics for both prey and predators (not shown) and 
huenwgeremsdependentprmobibity  saforiethe* epredator] ? (typere4) i 
Altogether, about 100 models were simulated, including some 
replications. Aussummary@oficesultssise@shownemnmfabie dad. 


A study involving refuge(s) for the prey was also performed. 
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5.1.3 Results 

Most of the simulation runs were carried out for 100 
time steps unless the predator or prey went extinct. In 
these cases, results were recorded at the time of 
extinction. 

Reproducing Murdie's experiments could be done easily 
using ECOSIM except for the satiation level which had to be 
set at 20 instead of 10. Prey varied between 10 and 30 
individuals and predators between 1 and 10. Some resuits 
were identical to Murdie's. For example, the number of prey 
killed was proportional to their initial number with both 
uniform and clumped distributions. Also, the introduction 
of a satiety level for the predator decreased the level of 
predation significantly. However, results indicate a higher 
level of predation when the prey were clumped, especially 
with a high prey density and when more than 1 predator was 
used. The fact that lower predation for clumped prey was 
observed in Murdie's experiments could be explained by the 
low density of ‘predator’ in his® studies. Tf a region 
contained only one or ‘two predators, and these predators 
were placed at random, it could take a long time before they 
reached the prey area with a random walk search. Especially 
when starvation was introduced in the model, a predator 


could easily die before reaching the clump. 
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Table 5.1 Results of simulation experiments 


on Murdie-like models. 
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When 5 or more predators were introduced, with no 
Satiety characteristics, the prey species went extinct in an 
average of 59 time steps. A typical Grunts ishown- “in 
Fig. 5.1. This surprisingly long period comes from the fact 
that if not enough prey were introduced, a predator could 
search for a long time before finding a prey. For some of 
the cuns, 20 successive failures to find a prey led to the 
death of the predator. In all but one case, the predators 
went extinct before the prey after an average of about 38 
time steps. The same results were obtained with the 
inteoduction of mobility for the prey (distributed. at 
random). 

When the mobility of the predator was made dependent on 
its satiety a slightly lower rate of predation was observed, 
especially for non-clumped prey, the predator staying for 
longer periods in "empty" areas. 

Interesting results were obtained by introducing 
breeding, first for the predator only and then for both the 
prey and the predator. In the first case, the predator 
species went extinct, after an increase to about 30 
individuals (from 10 initially), despite the breeding, 
mostly because of an insufficient number of prey (see 
DUC eos oe eto llce 1 tive second case, coexistence was achieved 
for at least 100 time steps, with both predator and prey in 
sufficient number. Values for the breeding characteristics 
were identical for the prey and the predator: a breeding 


interval of 5, a breeding probability of 1.0, and an 
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Fig. 5.1 Simple model without satiety for prey 
(species: l-prey, 2-predator) 
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Of ESpr ung distribution Of MCT litine With isan Initial 
environment of 10 prey and 1 predator, . 27 prey and 76 
predators existed after 100 time steps. It is likely that 
the predators would have taken over and annihilated the prey 
and then starved, but increasing the breeding values for the 
prey could lead to coexistence (see Fig 5.3). 

With the introduction of refuges for the prey, chances 
of coexistence increased. A system of type (3), with 20 
prey and 5 predators was run for 150 time steps without the 
extinction of the prey despite the fact that no breeding was 
introduced (see Fig 5.4). In this case, the refuge was made 
of 30 cells in the middie of the grid. When 5 small refuges 
of 6 cells each were introduced, the prey died faster than 
with one refuge, mainly due to the fact that a refuge would 


usually have a predator nearby. 


5.2 SYSTEM REQUIREMENTS AND PERFORMANCE 

ECOSIM iS written in ALGOLW, a revised version of a 
language designed at Stanford University. The ALGOLW 
compiler is available under the Michigan Terminal System 
(MTS) at the University of Alberta. MTSs) 1S, a. virtual 
memory, multi-task system developed at the University of 
Michigan to support both on-line and batch jobs. It runs on 
an Amdhal 470V/6 at the University of Alberta. A paging 
system is used for memory management. Other features of MTS 
include the supervisor which handles scheduling and resource 


allocation, Device Support Routines, and a spooling 
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Subsystem (HASP) which handles I/O operations. 

ECOSIM uses 90 4096-byte pages of memory during 
compilation, and the object file can be stored ina file of 
3224 lines. It takes approximately 8 seconds to compile the 
program, depending on the computer load. 

Parts of the program are written in FORTRAN and 
ASSEMBLER. Most of the FORTRAN routines are system routines 
used to interface with MTS. The random number generator was 
written in ASSEMBLER using a version created at the 
University of British Columbia. 

Thedresponsei setinmes vole rECOSIMe iss fastpradespite othe 
extensive dialogue designed for the novice user. Except 
when a large number of species is involved, or if the output 
interval is quite large, ECOSIM responds virtually 
instantaneously. 

When ECOSIM is run, 10 pages of memory are allocated by 
default. If the user is planning to create more than 800 
individuals or 1000 barrier cells, he should increase the 
number of pages to be allocated. 

It takes on the order of one tenth millisecond of CPU 
time to create a barrier cell and one millisecond to create 
an individual of a species. To process 1 mobile individual 
Fovcid00.timeststepsn (takeéess lapproximatelyo iat sseconds: A 
complete system of 10 mobile breeding prey individuals and 
10 mobile breeding predators with satiety characteristics is 
processed in about 1/4 CPU second per time step. More time 


is reguired when involved output is produced. 
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ECOSIM iS somewhat dependent on the MTS system for 
certain I/O routines and would require minor modifications 
to be made portable, but could be adapted to many large 


computer systems having ALGOLW compilers. 


5.3 CONCLUSION 

ECOSIM is designed to provide a useful and easy tool 
for ecological modeling. An involved dialogue provides the 
interactive facility needed for the model definition and the 
monitoring of the simulation. Realistic environments can be 
Simulated with the many available parameters, and even 
unrealistic studies can be done to discover the effect of 
parameter variations. 

ECOSIM could be improved in several ways. FLpSt sy) “ene 
environment is uniform. Each cell is identical except for 
the inclusion of individuals and barriers. No "physical" 
variables characterize the environment, like the moisture of 
the soil or the temperature. The individuals, as they exist 
in the system, are not distinguished by sex. Breeding 
specifications for a given species apply to all individuals. 
Thus all individuals are "female". Mating is not included 
in the  systen. Breeding is done at regular intervals 
whether or not an individual is alone, No age structure is 
included. An individual will not die or reproduce according 
to its age but according to the global function, predation 
and breeding variables. A function could be introduced and 


checked so that reproduction could be made dependent on age 
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and so that individuals would die accordingly. 

Finally, predators do not seek prey. In the system 
design, it was assumed that a cell is the physical immediate 
environment of an individual, as far as it can see and 
predate in the time step. 

The use of computers for simulation in community 
ecology is an enterprise still in its infancy. It may be 
hoped that availability of general purpose simulation 
systems like ECOSIM will advance both the study of ecology 


and thebart of ssinulatizons 
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APPENDIX 


USER'S MANUAL 


A.1 INTRODUCTION 

Before looking at the conversation between ECOSIM and 
the user, the procedure to use ECOSIM and the general 
constraints of the system are given. ECOSIM is easy to use 
and error recovery is simple. Most Lsofeithecetime;) the 
question is repeated until a valid answer is given. If the 
Situation is ambiguous, the error is pointed out so that the 


useGricanmecorrect it. 


A.2 GENERAL NOTES 

In the following, computer entries are in capital 
letters while user's entries are in lower case letters. The 
procedure to run ECOSIM is as follows: 


Signon i.d. 

PASSWORD? 

password 

PERE Tear. 

**LAST SIGNON ... 

GSE He Sista teres 

$rcun Mimiz:ecosim t=(number of seconds desired) 


te th FE HE ey He HE 


From then on, the user supplies answers to questions. 
Restrictions about the answers are as follows: 
1. No. answerseshould excecdoyZ20rechardcters: Additional 
characters will be ignored. 
2. No more than two numbers can be given at one time. 
Additional numbers will be ignored. 


3. Integers must be between -2147483648 and 2147483647, 
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Tals 


Absolute values of reals must be between 5x10-79 and 
7X107S, 

When a question ends with a question mark, the user can 
answer by YES, Y, NO,;*or N;: 

Some answers must be integers (like coordinates). The 
question will be repeated if a real is given. 

When an answer should be a real, an integer is 
conSidered valid (the decimal point is not necessary). 
If a pair of numbers is requested, they must. be 
separated by a blank. These numbers will always’ be 
integers. 

A temporary file is created at the beginning of the run, 
and all the information about the state of the model is 
stored. At the end of the run, the user has the choice 
of making the file permanent and/or printing the result 
at the batch printer. 

If plotting is requested, the user has to execute an 
extra step after the run is complete: 

$RUN *CALCOMPQ PAR=FILE=-PLOT 

Thestusen! idwsththeneiniormed*sofithe cost of@therplot and 
asked tfiehesstilliwishesetovget itis Tpeyes; tear earecei pt 
number is issued and the plot can be picked up later. 

If the user is planning on creating a big environment, 
he should specify PAR=SIZE=##P on the RUN command, where 
## is the number of pages required. The default option 


is 10 for about 800 individuals. 
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A.3 CONVERSATION 


In what follows, YES means YES or Y, NO means NO or N 


and REAL means a real number with or without a decimal 


point. 


EX PLALN 2o ese cites cd smlels ote ote oe ae oe oes BES OR NO 


If the answer is YES, the user is given a resume of the 
features of the systen. 


SRE LE evelelatsve ie cia ae alstetete eee olereatote cieicte ciate cain ee OR NO 


Tieetheotanswer \isPYES ASthe wsertwill) beaaskedvasserites 
of questions about the model he wants to create. Ef 
the answer is NO, a filename is reguested to retrieve 
an existing model. 


DESCRIBE ENVE RONMENT feats o oleletolexetelele e's ateteetelerele sete ctelee Volo OR NO 


EDIT, 


If the answer is YES, the user is given a resume of the 
Characteristics of each species and barrier and each 
individual of the different species and barriers at the 
terminal. 


SAVE, (G45 OR Ril Nive steve eletete sles e elclele efaleeslallir Sy OF OR R 


If the answer is E (or EDIT), a series of questions is 
generated to execute the required changes. 

If the answer is S (or SAVE), the user is asked for the 
name of a file which is created, and the current state 
of the model is stored for further use. The different 
variables are stored in a compact form not readily 
readable from the file. The file name is given the 
prefix OfSUSERS@9S0 thatethe) fuses tcaneizecognize, the 
filesMnet rcocatediduning theirun. | bForprexanple; cif the 
user creates a file name WORLD1 the real file name will 
be USER.WORLD1 but if the file name is requested (for 
retrieving information),,.the, modeler, uses, the namerhe 
gave to the file ise. WORLD1. 

If the answer is Q (or QUIT), the user is asked if he 
wants the temporary file containing the information 
about the run to be made permanent. Tt eyes) tea rifivde 
name is requested and saving is attempted. The user 
gets another chance if thesfile exists already or ‘the 
name given is illegal but not if there is not enough 
disk space to save the file. 

If the answer is R {or RUN), the simulation is 
initiated with the current state of the model (newly 
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created or retrieved from a file). 
EDIT, foAVE,ROUITS aOR CONTINUE? . cea seeks eteleteremiate co, MOF OP,OR RC 


Same as the previous question for E, S, and 0. If the 
answer Bris Cy athe Gsitnulation iseresumed@ where it was 
halted with all parameters preserved. 


RUN COUPLERS. SRE SULIS. DESURMD te seiela este less cisle cto ates ORONO 


Lfiethehanswer f{isteYEs, Gakabist soft allenspeciesiiand 
barriers with their characteristics and the 
Characteristics of each individual is given. 


SPEECHES WES. aye cael ete ere seal Grete cr aa atel Slate clave clas on ot wee emer OR NO 


If the answer is YES, species number N will be created 
(questions will be asked about its characteristics). 
If the answer is NO, questions about the barrier({s) to 
be created will be initiated. 


DISPLAY SYMBOL Scteatetee eielercielele ol stelelate ate cle clele cielelcieretetenl CHARACTER 


Any character except % is valid to represent a species 
Or barrier. When the environment is displayed, each 
individual of the given species or barrier is 
represented with this symbol. 


PROCESSING Ey tay ote et avelel cai cvahe 2c, sia el atatalvelolovebets lst cliete es cheleta ctoters ms GE 


A positive integer is requested. Motile species with 
the lowest processing level move first. 


MOBUEIT W(0ed) = Saws « alelel colar wialala’'s 3) cl acalny ebetal aia et eval. di eresaiecalorarer ersten 


A number between 0 and 1 (inclusively). If mobility is 
0, the individuals of the species never move. ff 
mobility is 1, each individual will attempt to move at 
Sveryintime® Stepar Thamobiiitytisetfor example 0.6, the 
probability of each individual of the species ~to™ move 
at every time step is 0.8. 


COMPETITION FACTOR (0-1) sc asccncscccesscceess ccccsccccocccs REAL 


AUSnuMber® betweeny ONSanrd. 185 (inelusively). Thi sss 
requested only if the species is motile. A competition 
factorweofys sayerd0s6,4 meansiathat if an. individual 
attempts to move in a cell that another individual of 
the same species occupies, he will "die from 
competition" with a probability of 0.6. 
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UNIFORM RANDOM WAL Kid Stetarer at stacsta cccra lela ere eleleler es A a a aac ES OR NO 


This is requested only if the species is motile. If 
ChesanswensistVES pkal probability ofel77 Mis assigned to 
each of the seven variables representing the six cells 
around and the cell occupied by the individual. 


6 WALK. PROB. (CLOCKWISE FROM UPPER RIGHT)....6 REAL NUMBERS 


Dhistis requested fonly it thetspécies istmotile: Each 
real number (between 0 and 1) must be on a new line. 
Diether Given. probabilities: are-a,), c;.0, se, and st. 
they @@areyhassigned,,clockwise,y from, -upper,. rights (see 
Farge wave Wi) The center cell (g) is assigned a 
probability of 1- (atb+tctdtetf). 
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Figs A.1 Assigned probabilitie 


In 
Tm) 
|O 
In 
| 
|@ 
IS 
{Qu 
10 
=] 


walk. 


If a given number is invalid, the user is requested to 
reenter it. If the total of the six probabilities 
given exceeds 1.0, the user is asked to reenter all six 
probabilities. 


PREDAMO RG telele ec eletetelelecie cteloloclelece ctetiele elelele erarelelaterc elevete statlilne OR NO 


Pee 


If the answer is YES, the user will then enter the 
predator characteristics of the species.) if the answer 
is NO, the given species will not have predator, 
breeding and hunger characteristics. 


SD BiG Sis) Scr evcren ee ccatiete a slelorelel cues) claret oete-so ohn deny Gib ty OR NOTHING 


A positive number or null =4Mline. If an integer is 
Given, a table 15 "created containing, the = predator 
Number and the prey number. When all prey species have 
been Specified for (thes givens predator, che) user 
continues by entering a null line. 


2 ; ra 


a. 7 ao 


ov ri PAGe cats Sede 


aparee 
2 «pl “as | 

ot her 2 ane 
ailae re Tes Sarivers 


aasanten die? seni A 6 si 


oe) aati sna ed 


anny sdgie 


(py 


heeoiensn 21 ee 
~ ae 

-_ a 

+-€ Per —e 

<e& 3#¢.9 2 ; . ao aa 

=—<t08 S-- j 7 

< ><a f > 7 

~=-< & p= - 7 Fk - 

_—— an 

a 


vibe webces tod cestei odessa tanptnas ras “on : 


of hersempe 2 qeeu add teloeet at sedge aan2e » ABO? OO 
wel?iiidesosq =le -aa to Intes edt TE oo eye Oe , he 


via Iie wtwees of Baves ep amen’ odd ae oe govip | 
Sidedomqy 
Oe 


of ou ANVrk » acces ede cde be awabe) ceekeeacenabedaceat(Gneneae - Re 


at 

ed? vetam fed? J Ske egy | — 
wmvann el? ei i : 7 - 
weicbasy «ved ton - 

apd n : 


sheet! os "heats 


83 


PROBABIEPTY OF TL LN Gite avahete sea etalaietotatet tn a cl ode ations Be 
A number between 0 and 1 (inclusively). The number 
denotes the probability that the predator will kill the 
prey if he encounters one. 


MAX KILLS / UNIT TVIsM Ehive eraveete tetatetars sie ete ele alate cle ate eteiee aie an LL GES 


A positive integer which indicates the maximum number 
of prey thatitheepredator facangaks!1t ad trcingtinene. tine 
step. 


BREEDING TEMES 6.5 elalelelelelolelecieleleleleleleie islets elelaleteralaislalsterstenlN ing ER 


A positive integer giving the interval (in time units) 
between breeding events for this particular species. 
Theteinitial Gage tofiethesiindividualstareidistributed 
uniformily over the age interval. 


BREEDING PROBA Biel ei Ve spapecstateeetelsteteteretoretete eeeeoee#esese AD OO BD OO 05 oD eb 


A number between 0 and 1 indicating the probability of 
a breeding event when the breeding event time is 
reached. 


OFFSPRING MiNi ebats (stowstel sfalereleralels atstelstelatets a (arsiatele sis stetotelereetete ere nll: 


A real number representing the mean of a Normal 
distribution used to determine the number of offspring 
to be created. 


OFFSPRING WARD ANCE Sete avalee cere chelcletetavelete ce cle cle tele fete: steleloteto snus 


A real number giving the variance of the Normal 
distribution used to determine the number of offspring 
to be created. 


SA ay. S GAtoe UOINO aA Reo OOO DOOD OO HO OOS OOo OO Ob OO Ooo 6.66 INE 


A real number indicating the satiation level of the 
species. mf ithe wsatretyrelevel oftan iindiuvidualsis 
equal to thetsatietyiscalestop, hesdeestnotteatheven eit 
he encounters prey. 


SA Ns SGARE BOLE OM sivarerer claceterevorels etebekereteretc atatetavarstarstateleteredoter lucas 


A real mumber less than the satiety scale top 
representing the starvation level of the species. If 
the satiety level of an individual is less than or 
equal to the satiety scale bottom, the individual dies. 
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+ BONS & Po Nal s aval a alain! 9 a%e cisva a7 al'cle loko’ Sneteo'a sl chelaie se ak Wate eae 


A real number which is added to the hunger level each 
time an individual of the given predator species eats a 
prey. The hunger level does not exceed the hunger 
Scale top! forthe species. 


ay ENCRER ENG: Qalate ste « Scie a cletcls sierein tee eae tine erent. One A OREAS 


A real number which is subtracted to the hunger level 
each time an individual of the given predator species 
Cats a prey. 


S Averys JEM MERRIE NIE MEVAN SMivovsracaieo%s(e alata ave taretalere co ctetetatstanel ete sate akeretn ad: 


A real number representing the mean of a Normal 
distribution used to determine the initial hunger level 
o9f each individual of the given species. 


SY AME IIE VE i Nae eA, VA ROSA NIGH Shovetele clotetee clele stele croreteteterete clelelers ebsrer cal, 


A real positive number giving the variance of the 
Normal distribution used to determine the initial 
hunger level of each individual of the given species. 


SATLIECYRDERENDENT MOB UUETY Goes cistels se tie soe ese ee cece st us OR TNO 


TEAthe,sanswer.is,YES,,theseprobability .of, moving for} fan 
individual at a given time step is determined not only 
by@his!mobility facter buttaksotbyy his vhungere levels 
The probability value £(x) will be on the curve passing 
through (0, HUNTOP), (MOB.COEFF., (HUNTOP+HUNBOT) /2), 
and (1, HUNBOT) where x is the hunger level of the 
individual, HUNTOP and HUNBOT are the hunger top and 
bottom, and MOB.COEFF is the mobility coefficient of 
the individual. 


eX OU) Sa ON Sa erie atela auals slelelele ea slale ulais! e Sontere ba scala emia eo are ee OR NO 


If the answer is YES, the user is asked which species 


arerexcluded :(ises which! (species) rindividuals )cannot 
move in a cell occupied by an individual of the given 
species. 


SPECIES NUMBERS... 2a2ees JES ES. TSEC es NS es Mstane al NTSGER 


An integer greater than 1 specifying the species to be 
excluded by the active species. When no answer is 
given (null line) it is assumed that no other species 
is excluded. 
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NATE AL NUMBER Sees o)6)'s) ele 4.8.0) sels ehelai ciel alse ale ists 5 s/avaleletele eters cuen DhGak 


A positive integer representing the initial number of 
individuals to be created. 


RANDOM Daa BD ie) Noon ot ctaet a area aia create onan a a aoe ES OR NO 


If the answer is YES, individuals or barrier cells are 
distributed according to a random distribution (chosen 
by thevuser after); af ‘the answer 1s) NO, - specific 
coordinates will be requested. 


U, N, Ey OR PD Sfalskateda elareterete ate le a aterets o cinletstalaraetertre ate) CHARACTER 


The answer should be the character U, N, EF, or P. igs 
the answer is U, the individuals or barrier cells are 
distributed according to a Uniform distribution in X 
and a Uniform distribution in  Y (parameters are 
requested subsequently). If the answer is N_ the 
distribution used is Normal, if E, Exponential, and if 
P, Poisson. 


X(Y) NIMPEN AUN voreletonere lek ouelel erelevcichelelelsrersteteleveroretebete el aterccekelenerehotereneieren cur nny 


A real number representing the mean of the Normal 
distribution of the X(Y) coordinates. 


RCV NVAR? oa 2 oso 00 0 Blolelolatelelotetetelelcleteloletelcieneiataletercteleleleteteleteterenater cnt. 


A positive real number representing the variance of the 
Normal distribution of the X(¥Y) coordinates. 


RAV) NLIMENE (SUP) 20. saa ss slsle © o'n'a610 0 (aia ee 4/4 01010 6 «08s © o/s 60 e slain 


A real number indicating the limit for rejecting values 
generated according to the Normal distribution. 


X(Y) MIN ca csw cna ssa cesses enesssesssesse Relea acl etele diel ete si clotsreskr Ads 


A real number indicating the minimum value for the X(Y) 
coordinates of the individuals distributed uniformly. 


X(Y) MAX 2 oven e cae ccccccccnenssessscesesesesssscesessccccs REAL 
A real number > X(Y)MIN indicating the maximum value 
for the X(Y) coordinates of the individuals distributed 
uniformily. 


Kes AM Bisteneteletalelslere easoeaaeaeneananvvaensvaeeveeseewnavan ea as Sate Sb etels Sales ste ROAD 


A real number giving the value of lambda in the 
exponential random distribution. 
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XEORIG 


shele'e eielele O10) 0), 056: 4/0 © O46) 6 6) 1 61610) 6 6) «1 0) 6) 6 aa 014) 4.4 4 ale of cles elelareteta enna 


A reai number indicating the origin of the exponential 
curve. 


tL UM Detect tateberet stele! o°e o shad ctete a raievatatetata cco cleat, es eee. ea AL 
A real number indicating the limit for rejecting values 
generated according to the exponential distribution. 


The following variables are similar to the previous 
Ones, but apply to different distributions: 


YELAMB, YEORIG, LELLALES XPLAMB, XPORIG, APLIOINE, 
YPLEMSUP, “YPEAMB;, YPORIG, VPLIMINE, aL 2Linsvue 
N X-Y ACOORDINATEOSP APRS si. 2 cle s slalsre c'slns wee see allege. NTEGODREPATRS 


N integer pairs should be entered, one pair per line, 
undicating the "positions of “the grid in which the 
individuals are to be placed. Or, the user can enter 
S$CONTINUE WITH FILE RETURN, where FILE contains the 
given number of coordinates, one pair per line, each 
humber separated by a blank. 


PRENTPINLITGAL DISTREBULGON Ss o cle ses clu ss 96 sie celteses sie o LLO pORAND 
TE the anSwerobsaSeeves, oathe  edistributbionl eof oathe 
individualsoist#listedg@in the fonmpoindgae tablethgi ving; 
fOLeReaAcChe individual, jits, poszoront(s, Yjytitstage, its 
hunger level, and the number of prey killed. 


GLOBAL MODINE MUON Sicietetets shel ststsletele slotaletstelsnetototolstetatsteleie Lo OR NO 


If the answer is YES, a function is then requested that 
is applied to the population of the given species at 


every time step. If the applied function implies a 
decrease in the population, individuals are deleted at 
random. If the result is an increase, individuals are 


cneatedgtaccordingumtoy the ¥einitizalfidistribution #£er 
position and hunger level. 


PUNCTOLON S20 Aa. 2 he Se Sake Sole chalc b aleleala sl shtele gUE TO 80 CHARACTERS 


Ae’stting | ‘representing. .asfunction.of,xX.wheresx9isrethe 
present population (number of individuals). For 
example X-2 Or 2+ (X**2). The string should not contain 
intervening blanks as a blank indicates the end of the 


string. 
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BARRIER N? 


s( Se (610) @lerele @ 14) 06 ofa) 6h eee 1066 016-616 6 s 416 0 clele eee ato OR NO 


If the answer is YES, barrier number N is created 


(questions aves tasked merorercharacteristics)- If the 
answer is NO, »vereationlfofe the imoded iiss. considered 
complete. 


SECT E Sitveesta cheleldlelelc.ea o eramioaie o malt awe bh eae wen LUTEGCER OR NOTHING 


A positive number 0 or nothing (simply press ENTER). 
If an integer is given, a table is created containing 
the barrier number and the species number. When all 
Species numbers have been specified for the given 
barrier, the user enters a null line. 


CROSSING PRO BARI LEB VGs «0 a che wo Siolsll o e.eie epaiete « inlateld Siele ws oats Se REL 


A number between O and 1. The number denotes the 
probability that the species individual will cross the 
barrier if he encounters one, 


NUMBER OF GEG Si ais cise o abate eo clelcle sloteFerc calc sla steleberelele cl ceters ONC EGER 


A positive integer representing the initial number of 
individual barriers to be created. 


SPECIE € (Sy OR CORNERS (Cee Selsletatae cclalela eck weleie cle caren, OR ‘é 


The user has the choice of specifying the barrier cells 
one by one or by giving the coordinate of the four 
Sides of a rectangle (left, right, top, bottom). 


LEPEA (RIGHT PRTOR M BOTTOM) Ss ces ec cnc s cence en eece sees gee UNULGER 


The LSEERE(LVGnt, (top, .pbotton))) «coordinatel sorer tic 
rectangle formimg the barrier has to be specified. 


PAcTeRN iM boner a a Nonlale atc o ete etarelerstetalsl anata siete wlelel silitts TO 20 CHARACTERS 


A file name has to be given. If the user requests an 
environment previously saved on file, the file must 
exist. If the file name is requested to store an 
environment or the results of a run, the file must not 
exist and is created by the system. 


TRY ANOTHER it Ties ate shavelarerekel ete MiGhala erate ee siwiaiele « scetelavetest cis OR NO 


This question comes if an attempt is made to create a 
file’ which. .already.exists.or,.to get«a.file.which do¢es 
not exist. If the answer is YES, a second attempt is 


made. 
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PU TPE ORR R ee ce ECCLES, ORs ENDi< cee welne wolg o ca saaeacH ens , OCR RE 
If the answer ispeByeeth@. «USGL sii Se psaskeds Tto enter 
information about the barrier he wants to edit. If the 


answer is S, species editing is initialized. nt OR ye 
answer is E, editing is terminated. 


EOD eve eG eee NOD TTY Aa SPHCIE St oets asec sere s)sate oe Onn 


Depending whether the answer is A, D, or M, questions 
are asked to create, destroy, or modify a species. 


WHICH SF ee PIG stata, atal ae ata oi abe arab ct crate aie oa e cree 2 ou a dlale antec EN DEGER 
If the answer to the previous question is D or M, the 
user is asked to specify the species (positive 


integer). 


BODIE Ys 
eo ous PLUAY SViICBO Is ive aretetelercteke:s eceelele atele.6 a wie ee eletece a tater eS OR NO 


If the answer is YES, the user is prompted for a new 

display symbol. 

Fora¢gthet thiol bowing pt Lisetheteusericanswers YES, hevas 

prompted for the change. 
sie SEROCESSING LEVEL is deste a ssinlels aclencccsassecevases tuo OK NO 
os es HOBDULITY RFACBOR Betas cca dg aiclsis ale slaves cscs scieessisisces tLhuaeOkR NG 
oie ew COMPETIVIONGEFACTORG os. te sidaacdlieosiceescaccenceseccse lino OR NO 
en sMALKINGTPROBABLITDIES? see sccacasdacsaee tases suseces Lon -ORPNO 
one PREDATOR CHARACTERISTICS Seis eke a civesis ceisicea sas LOO PORNO 
sists MAXMRGELSOY UND TOLUENE Ts 6 cle a oc'c'e a sleeve Scte eialele Se estar HO AGK. NO 
52s PROBABLDDEY POPPKILLENG OO. Sess Seca ces caesecsssettis OR NO 
SUOBREEBDING, INTERVAL oes ce sececiwecocedcaeenegeden cl Lo -ORRNG 
oe e BREEDING PROBABILITY? es csesaveoccrscssceaccssesee LES OR NO 
A, DOF FSBDRENGMMEAN? oo Stele « Sale sacle bide s Mats aS Sle omele ne ak oo BORENG 
PONOME SPRING VAR LANCE? cs occ ee cee sald cincl os esicnrenedee ens Ole NG 
Pee eun TMT YEPOR v4 esis aieialscletslelscisiejets else's else 9's 6a 0 7ee'ee et Cen nO 


rare ALL le Belay BOMUO Ml os 6 clea ets ateretbiste eleleo «ies ete aiae ste) eiaterely i) OR NO 


«eeHUNGER MOBILITY FUNCTION? cae cevcccsccscccccceee YES OR NO 
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Serer NC MIEN TY se @ «16ers sielelasisisis.e.e.0 0ig.eipahe's wie 4 cake roy eomO RONG 
sees ENGREMEN D2. 5s. Sralciste statetela c gavclele ses HAM. Ona PT AVES SORMNO 
ee HAG LUSLONS CHARACTERISTICS 740 eset cles wise seen ce see VES OR-NOG 
maroLO GAUUBUN CRTON 2 its. c ate gie'a es é/elasd odie «ica Celene oa cE OnOnn NG 
sie SU REY SS PECT ES tttevel sini cisc cc aidis oo ble o ahile's oa Tats « Aelnaeay ESUORENG 
cate oN DEVIDUAL S Actict sc stolec ac dhatelevels’ ol eisiielelsiciele ere ciste o ae ESMORENO 
see ol POTES/BARRT ER@CHARACTERTSTICS 224. Ge.03 700. ee SLES ORD NO 


aa COORDMS (C)rPESATIECYI“S)§, GORGBREEDINGOUS)>......Cp Sp OR 8 


Depending on the answer, the coordinates of an 
individual are altered, or his satiety level or his 
age. 


DOES SPECIES NB BECOME PA “PREDATONG@s a ste cies s sie eeisletetety (oy PORENO 


This information is requested if the user asked to 
Change the predator characteristics and the given 
Species waS not a predator before. 


DOESMOPECTES NyREMAINe Ay PREDATORS vise a ae wee se ee Lb BORG NG 


This information is requested if the user asked _ to 
change the predator characteristics and the given 
Species was already a predator. 


CHANGE, DELETE OR (NDI vane ele el avelelalelalaleleleleletetelelele) als) eelsra. Cig, Dy OR A 


If the answer is C, changes to existing prey species 
are initiated. If the answer is D, species N will no 
longer be a predator for the species specified. If the 
answer is A, a new record is created since species N 
will now prey on a new species. 


ADD, DELETE, OR MODIFY INDIVIDUALS: ....2.2eceeceeeeee A, D, OR M 


Depending whether the answer is A, D, or M, questions 
will be asked to create, delete or modify individuals. 
"Individuals" here is taken as an individual of a 
Species@ordan individual tbarrier cell. 
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HOW PUArty Lpals ale rah or 5 4's 0: a. 6a a ahaie etal evalahor et ete, aon aias We atk ow aot eT Pore 
A poSitive integer specifying how many individuals are 


to be created, deleted or modified (depending on the 
previous answer). 


CoV ameCOURD UN ATE DAT Rut sicsssmajete ter eietasstelele cela ata alae LN LEG E Rab A 
A coordinate pair has to be given indicating the 
position to add, delete, or modify the individual. If 
the user does not remember the existing individuals' 
POSttions7; §ehe. can presssEND OF EFIUESTOA. listing of the 
position of each individual will be given so that when 
the user is reasked the question "...INDIVIDUALS?", he 
knows where to delete or modify then. 

NEW COORDEN A DES 2s Melee co slelc a abe ste 6c0e Siew delet «cles DUN DE GER PAIR 


This new coordinate pair becomes the new position of a 
Umodrvived! individual. 


ADD, DELETE, GORWMNODIFY WASCBARREER Sees Selsces se oseslc shy Dy COnEM 


Depending whether the answer is A, D, or M, questions 
are asked to create, delete or modify a barrier. 


WHICH I URMRUE RO EA OO OO tur ode wate wvatavele shetets ete. ec etote scetetee ote Nee Gink 
If the answer to the previous question was D or M, the 
user is asked to specify the barrier (positive 
integer). 


UNC WHOECH PIME STEP! ss csc ses esses s celsicgscccese vice INL UGER 


An integer number specifying until which time step» the 
run is to be carried. 


OUTPUT HAN Ep RaVeAtistetete larolelelel elelalelelelere ee Slolelerclelelstsienereieveteteimere tl Nila Grune 


An integer number indicating the output interval for 
which the environment can be displayed. 


STOPPING CONDITIONS: SPECIES. cccescsssccccccccscccce NIEGER 


An positive integer specifying for which species do 
stopping conditions exist. 
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RU AAG Pat oi oh ai el'af aval sf o ap bias of oh a¥aiak abo of onchaiev cla? aren AEC eS PATR 


An integer pair giving the bounds on the population of 
the "given “species. ) Ifthe number of individuals’ of 
this species becomes less than MIN or more than MAX, 
the Simulation run is halted. 


RESULT Ss LEE 3B ea MADE PERMANENT? Joa oes ss5 oo aateee lolho Oriana 


If the answer is YES, the temporary file which contains 
all the results of the simulation run is renamed to be 
made permanent. A filename is requested. 


PRINT RESULT PLL Pe istekstetetetelalc a's a otetale ter erethietcla tie ots cre aS or NO 


Tf the answer is YES, the file contdining the result of 
the Simulation run iS printed at a batch printer and a 
receipt number is issued. 


PlOpl ie teeceia alaleiece «eas blela’ alee 616 (os ale-ats os) 6 sts ele 6 ete e eeteteee he (Oy & NO 


If the answer is YES, species numbers are requested, 
For each species, population versus time step is 
plotted on the Calcomp plotter. Ali species are 
plotted on the same graph. The user can separate the 


plots on different graphs by answering YES to 
successive "“PLOT?"s, At the end of the run, the user 
Should 


SRUN *CALCOMPQ PAR=FILE=-PLOT 
if he wishes his flots to be executed. 


MAXIMUM S IGA) OF GReUD oc retoteletelalslonetelebeleletcrelslcletetaictcleleteralecereio cl Nari 


GRID 


An integer between 1 and 16 indicating tha xX and Y 
dimension of the grid (number of cells horizontally and 
vertically). A size of 14, for example, will create 
14x14 windows. 


(G} OR NO GRID (Nis erase ale clei e as 06 01 ele 1s) s 0101416 5 6ia'e 4.46761 or N 


Peet ieee answeCLe !SanG,y, atau Ulid 1s. sala Ri alountdsEene 
individuals when the state of the model is displayed. 
Tf the answer is MW, only the contour of the given state 
of the model is drawn. 


DISPLAY OUT Dill 2s os ateielefa slates ave esis) sisie clels/e ere elereersiete see) io or NO 


At each time where the output interval is reached, the 
user has the option of cancelling the display of the 
state of the model on the terminal (the state of the 
model is still recorded in the temporary file). 
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REGION: 


OC OVCHO ROO OM vO OMe yO iC Cu Or. Onin Oo olen OOO Ob o Cee np LORMEIn eine 


The state of the model environment is divided into up 
to 16x16 cells. The region numbers run horizontally. 
The user sees only the desired regions on the screen. 


tunel oC Vatatarsielaletereletelers oiatcleisreteta rene ee eee RS OR NO 


If the answer is YES, statistics are given about every 


species and every barrier for the present state of the 
model. 
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